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We study the fate of the Ising model and its universal properties when driven by a rapid periodic 
drive and weakly coupled to a bath at equilibrium. The far-from-equilibrium steady-state regime 
of the system is accessed by means of a Floquet mean-field approach. We show that, depending 
on the details of the bath, the drive can strongly renormalize the critical temperature to higher 
temperatures, modify the critical exponents, or even change the nature of the phase transition from 
second to first order after the emergence of a tricritical point. Moreover, by judiciously selecting the 
frequency of the field and by engineering the spectrum of the bath, one can drive a ferromagnetic 
Hamiltonian to an antiferromagnetically ordered phase and vice-versa. 

PACS numbers: 


The Ising model is undoubtedly the most studied 
model of statistical mechanics. Besides its equilibrium 
properties, its coarsening dynamics following a tempera¬ 
ture quench from the paramagnetic to the ordered phase 
is also quite well understood mm, even in the presence 
of weak disorder GHg. Taking into account the dissi¬ 
pative mechanisms due to the inevitable coupling of the 
spin system to an environment has been successful in 
the description of important many-body phenomena such 
as the decay of metastable phases IZHE], hysteretic re¬ 
sponses mm and magnetization switching in mesoscale 
ferromagnets OIlCEEj- As it is becoming clear these days 
that driven-dissipative physics, i.e. the balancing of non¬ 
equilibrium conditions and dissipative mechanisms, is a 
promising route to achieve a new type of control over 
matter, a burning question arises: can the Ising model 
be driven to non-equilibrium steady states (NESS) with 
enhanced or even novel properties? 

This question has been approached in the context of 
slowly oscillating drives (magnetic fields or electrochemi¬ 
cal potentials) by means of Monte-Carlo simulations [T5}- 
HEII25H2ZI, mean-field treatment [T9H24] , or other analyt¬ 
ical techniques p25H5T] . One of the key results is the ex¬ 
istence of a so-called dynamical phase transition, where 
the cycle-averaged magnetization becomes non-zero in a 
singular fashion. This has recently been supported by 
experimental evidence in the dynamics of thin ferromag¬ 
netic films (32l . 

In this Letter, we focus on the Ising model driven by 
a rapidly oscillating magnetic held hcos(cot). We depart 
from the usual Floquet engineering approach to many- 
body phases, mostly directed towards cold-atomic sys¬ 
tems I .T i l , by including a dissipative mechanism, namely 
by weakly coupling the system to an external equilib¬ 
rium bath. The properties of the latter are kept generic 
in order to study the influence of its spectrum on the 
dynamics. Dissipation is the natural counterpart of driv¬ 
ing and our results are an unquestionable proof that the 
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Figure 1: (color online) Mean-field magnetization \ip\ as a function 
of the temperature T of the bath (super-Ohmic, s = 1) for different 
values of the drive h given in the key: h = 0 (equilibrium), h, ui, 
h oj and h uj. The critical temperature for h -C u, T c = 
(1 + s)T® q , is computed exactly in Eq. jl6| . In the temperature 
range T c eq <T <T C , \tp\ ~ (h/u) 2 / s u> / zj. (w = zj = T c eq ). 


presence of a bath can have far-reaching consequences. 

We access the non-equilibrium steady states by means 
of a Floquet mean-held approach. We derive the mean- 
held self-consistent equation for the magnetization and 
use it to derive the non-equilibrium phase diagram. 
Whenever analytical solutions are beyond reach, we com¬ 
plete the picture with numerical results. Our main results 
are to show how to combine drive (i.e. h and w) and 
dissipation (i.e. mostly the low-energy spectrum of the 
bath) to increase the critical temperature T c , to modify 
the critical exponent fix, as well as to change the or¬ 
der of the phase transition. Additionally, we show that 
the drive can, in the presence of carefully selected baths, 
convert a ferromagnetically ordered system to an antifer¬ 
romagnetic order, and vice versa. 

Model. The total Hamiltonian is composed of the 
system, the bath and the system-bath Hamiltonians, 









2 


= Hs ( t) + H b + Hsb 

with (we set h = k b 

= 1 ) 

Hs(t) 


- h cos (cot) 55 a i ) 

(la) 


(ij) 

i 


H b 

= 55 0Ja b i,a bi > a 

i,a 


(lb) 

Hsb 

= 55 ta a i (\a 

i,cx 

+ b i,a) ' 

(lc) 


The spin 1/2’s, represented at each site i of the bipartite 
lattice by the usual Pauli operators a*’ y ’ z , are interact¬ 
ing through a nearest-neighbor interaction J. h is the 
strength of the periodic drive with frequency ui = 2tt/t 
(we choose co > 0). Equilibrium conditions are recovered 
for h = 0 or co = 0. Notice that in the absence of an 
environment, the drive has trivial consequences on the 
dynamics of the Ising model. Indeed, as [af,H(t)] = 0 at 
all times, all the degrees of freedom are conserved quan¬ 
tities. The environment is composed of local baths ex¬ 
pressed in terms of a collection of non-interacting bosonic 
modes labelled by a, with energy co a and with creation 
and annihilation operators given by b\ a and Each 
is in equilibrium at temperature T = 1//3 and we as¬ 
sume it is a “good bath”, i. e. it has a very large number 
of degrees of freedom and it remains in thermal equilib¬ 
rium. Below, we replace '}2 a by f dep(e) where p(e) is 
the bath density of state. Without loss of generality, the 
chemical potential is set to 0 and p(e < 0) = 0. Hsb is 
responsible for thermal spin flips and the parameters t a 
control the strength of the spin-bath couplings. After in¬ 
tegrating out the bath degrees of freedom, the bath will 
enter the reduced problem via the hybridization func¬ 
tion u(t) = \t(e)\ 2 p(e). The low-energy behavior of the 
hybridization i/(e) ~ e 1+s characterizes whether the 

bath is Ohmic (s = 0), sub-Ohmic ( s < 0), or super- 
Ohrnic (s > 0). Note that we do not consider additional 
system-bath coupling terms such as af’ z (bi^ a + b\ a ) be¬ 
cause they do not induce any qualitative change in the 
non-equilibrium dynamics. 

Floquet mean-field description. The time-dependent 
mean-field Hamiltonian corresponding to H (t) in Eq. (jTJ) 
is the one of a single spin coupled to its local bath, and 
reads H(t) = Hs ( t ) + Hb + Hsb with 

Hs {t) = ~ zJp ( t ) a z — ha z cos (cot) , (2a) 

H b = 55 b'l b a , (2b) 

a 

HsB=^t a a x (b a + bl) . (2c) 

a 

Here, p(t) is the expectation value of cr z (t) which serves 
as the order parameter, and z is the coordination num¬ 
ber of the bipartite lattice. When the coupling the bath 
is weak (see the discussion below), the spin subsystem 
can be seen as quasi-isolated during many periods of the 


drive. There, the Floquet theorem states that the in¬ 
stantaneous eigenstates of the time-periodic Hamiltonian 
Hs(t) can be written in the form \if a (t)) = e _1 ' E “*|')/’ P (t)) 
where E a is a so-called Floquet quasi-energy and |^/> p (t)) 
is periodic: |V’a(t + T )) = |f/> P (t))- Owing to the fact that 
a z is a conserved quantity, we may choose our Floquet 
eigenstates to simultaneously diagonalize a z . Note that 
this also implies that tp(t) is a constant (at least between 
two events induced by the weakly-coupled bath). Alto¬ 
gether, the instantaneous eigenstates of Hg(t) are simply 
given by 

|t(i)) = e + i [ zJ W+S sin( “ t) ] |t) =e -i|f p (f)) , (3) 

\l(t)) =e - i [ zJ v t += sin (“ Jt )] | l)= e - i ^ t || p (f)) , (4) 

from which one identifies the Floquet quasi-energies and 
the periodic states, reading 

e t = -z J<p , |t P ( t )> = ]T Jn (h/co) e~ inu>t It) , (5) 

n 

e l = zJp, |t P (0> = 55 Jn(h/u)e +inu,t \l) , (6) 

n 

where J n are the Bessel functions of the first kind. 

Rates. The bath is inducing incoherent transitions 
between the eigenstates. The transition rate R^i from 
|t) to |t) can be obtained by means of a Floquet-Fermi 
golden rule [5]: 

Rfliv) = 27r 55 l A tII 2 9 ~ e l + > ( 7 ) 

m£Z 

with g (e) = ^(e)[l + «b (e)] + v(— e)riB(-e) where the 
Bose-Einstein distribution ras(e) = l/(e^ e — 1) and 

r T Af 

Aft = / - (t P {t) | a* |t P (t)) e imuit = J m (2 h/co). (8) 

Jo T 

A similar expression can be obtained for -Ri-|-((/?) with 
A™ = Aft". Note that the integration over the degrees 
of freedom of the bath also contributes to a small renor¬ 
malization of the spin Hamiltonian (so-called Lamb-shift) 
that we neglect. 

Steady-state population. We stress that the previous 
analysis is valid only in the case the bath is weakly cou¬ 
pled to the system, i.e. the rate at which it induces spin 
flips is much smaller than the frequency of the drive: 
Rfl, -C co. Under these conditions, ip(t) is indeed 
constant over many periods of the drive and a time- 
translational invariant non-equilibrium steady state can 
establish. Once it is reached, the probabilities of being 
in the | f) and | 4-) states are simply given by 

P ” ESS = ITfibiCf ‘” d ^ = 1 “ P "“ S ' (9) 
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Self-consistency condition. The probabilities in 
Eq. @ allow to compute the steady-state average 
magnetization as \ip\ = |P^ E8S — pN ESS |. Therefore, we 
obtain the self-consistency condition for the mean-field 
order parameter 


Here, the + sign corresponds to a ferromagnetic order 
while the — sign corresponds to an antiferromagnetic or¬ 
der. Making use of the expression for the rates given in 
Eq. 0 , we obtain 


Rj ~ fl-nM 

RM<p) + P t -iM ' 


( 10 ) 


Ri t(v) - Rniv) = 

Rtf(tp) + R n (ip) = 


2n | J 0 (2h/uj)\ 2 v (\2zJtp\) sgn (Jip) + 2n EE | J n (2h./oj)\ 2 bv (anu; + 2bzJp) , ( 11 a) 

n>0 a,b=± 

2n | J 0 (2h/u>)\ 2 v (\2zJp\) coth {j3\zJp\) 

+2-7T EE | J n (2h./u>)\ 2 v (anu) + 2bzJip) coth(/3(a?rw + 2bzJtp)/2) . (lib) 


n >0 a,b=± 


In case the ac drive is switched off, h = 0, one naturally 
recovers 

Ri i-M ^ Ru(v) 


±p = 


-> tanh fdzJp 


( 12 ) 


RaW) +Rn(v) h =° 

which is the familiar self-consistent condition for the Ising 
model in thermal equilibrium. In this case, it is well 
known that there is a second-order phase transition at 
the critical temperature T® q = z\J\, below which fer¬ 
romagnetic solutions are possible for J > 0 and anti¬ 
ferromagnetic ones for J < 0. 

Non-equilibrium steady-state phase diagram. The 


self-consistency equation (10 1 together with Eqs. (11a) 


and (lib) allow us to explore the complete mean-field 
phase diagram far from the equilibrium regime. Let 
us first investigate the fate of the well-known second- 
order phase transition in this out-of-equilibrium con¬ 
text. In order to access its locus in parameter space, 

Using 


we expand and solve Eq. (10) around ip = 0. 


the low-energy parametrization of the bath hybridization 
v(e) ~ 77 e 1+s , we obtain 

e— >0+ 


±p = 


Rj ~ Rniv) 
Ra(<p) + Rniv) 


= fdzJp 


K\2zJip\ s + A 
K\2zJtp\ s + B 


(13) 


where 


K = iq \J 0 (2h/uj)\ 2 , 

A = 2 | J n (2h/u>)\ 2 v'(nco) 


n> 0 


B(T) = PY, 


n> 0 


Jn\~ 

OJ 


v ( nco ) coth 


/ fdnui 


Besides the trivial solution tp — 0, the self-consistent 


mean-field equation (13) admits non-zero solutions 

1 f B(T) ±sgn( J) [A/B(T)} T® q - T] | /s (M) 


2 T r eq 


I< 


T T sgn( J) T c eq 
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Figure 2: Non-equilibrium phase diagram in the drive vs temper¬ 
ature plane for different values of u given in the key. (.s = 1). 


Equation (14) above is quite rich and its analysis below 


will tell us about 1 ) the critical temperature, 2 ) the na¬ 
ture of the ordered phase (and the stability of the non¬ 
trivial solutions), 3) the critical exponent, and 4) the 
nature of the phase transition. 


Note that ip in Eq. (14) must vanish continuously when 


crossing a second-order phase transition. For a bath with 
a sub-Ohmic low-energy behavior, — 1 < s < 0, this im¬ 
plies that the corresponding critical temperature, T c , is 
identical to the equilibrium case: T c = T® q . There¬ 
after, unless stated otherwise, we shall focus on baths 
with a super-Ohmic low-energy behavior, s > 0. In 
this case, the critical temperature is the non-trivial so¬ 
lution of T c = ±sgn(J) [A/B{T C )\ T® q . Before solving 
explicitly for T c , one can already remark that T c must be 
larger than T® q so that the numerator and denominator 
of Eq. (141 have the same sign for T® q < T < T c 


en¬ 


suring a well-defined non-zero magnetization solution in 
that temperature range. Let us now solve for T c by con- 
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sidering the case when h <C w, for which only the n = 1 
mode contributes significantly (because of the stronger 
power decay of the Bessel functions for larger n’s). In 
this case, the critical temperature T c is determined by 

tonh (2%)= ±sg,,(J) 2 w (15) 


Note that Eq. (15) has a finite solution only if the norm 


of the right-hand side is smaller than unity. 

Importantly, when i/'(u;) > 0, the type of order is dic¬ 
tated by the sign of J in the ordinary way: J > 0 for 
a ferromagnet, J < 0 for an anti-ferromagnet. How¬ 
ever, it is noteworthy that driving can turn a ferromagnet 
into an anti-ferromagnet and vice-versa when zA(u;) < 0. 
The choice of sign in Eq. (151 that yields a solution 


(phase transition) in this case is the opposite of the com¬ 
mon Ising model: here when J > 0, there is an anti¬ 
ferromagnetic solution, and when J < 0, there is a fer¬ 
romagnetic solution. 


Eq. (15) can be solved analytically when the right- 


hand side of the equation is much smaller than unity, 
v{uj)/\ v' {uj)\ <C T® q , yielding the critical temperature 


T c «T c eq 


(16) 


Eq. (16) transparently elucidates that by judiciously 


choosing the driving frequency or engineering the bath, 
or both, one can achieve a rather large critical temper¬ 
atures T c , much larger than the one for the undriven 
system, T® q . To exemplify this point, let us assume 
that the low-energy energy behavior of the hybridization 


'(e) 


,1+s 


(s > 0) holds up to the scale to. This yields 


T c ss (1 + s) T® q > T c eq . See also Fig. (1) where we plot¬ 
ted the magnetization as a function of the temperature 
for different drive strengths. In the temperature range 


T® q < T < T c , it can be seen from Eq. (14) that the 


drive is responsible for a finite magnetization on the or¬ 
der of \ip\ ~ ( h/uj) 2 / s uj/zJ. In Ref. [37], we show the 
stability of this non-trivial mean-field solution below T c . 
In Fig. (2), we summarized the non-equilibrium phase 
diagram in the temperature-drive plane by numerically 
solving for the critical temperatures in all the regimes of 


h and w. Beyond the super-Ohmic case, Eq. (15) suggests 


that one can engineering very high critical temperatures 
by using the edges of the bath spectrum to realize very 
large \ v'{ uj)\ or by embedding the spins in optical cavities 
with a finely tunable sharply peaked spectrum. 


Equation (14) also readily provides the mean-field crit¬ 
ical exponent for the order parameter as function of tem¬ 
perature, /?t = 1 /s, to be contrasted with the undriven 
case where the mean-field exponent is /3/ q = 1/2. This 
means that, even at the mean-field level, driving changes 
the nature the phase transition. 


Finally, Eq. (14) predicts a diverging magnetization 
at T = T® q . Although it was derived under the as¬ 
sumption that kp is small, this suggests that the original 



Figure 3: Non-equilibrium phase diagram in the temperature vs 
drive frequency plane for fixed h/ui = 0.1. The red circle indicates 
the location of the tricritical point separating a second-order line 
below from a first-order line above. (s = 1). 


self-consistency Eq. (10) may have non-trivial solutions 
kp / 0 which are not connected continuously to kp = 0 
and signaling the presence of a first-order phase transi¬ 
tion. For example, in the case of baths with a sub-Ohmic 
low-energy behavior (—1 < s < 0), the denominator of 
Eq. (10) given in Eq. (lib) has l/(kp — ip n ) divergences 
located at every kp n = nu>/2zJ for n = 1... [2zJ/lj\. 
In turn, this implies the presence of a collection of non¬ 
trivial solutions of the self-consistent Eq. (10) close to 
these kp^ s. For baths with a super-Ohmic low-energy 
behavior, the denominator Eq. (lib) is well-behaved and 
we investigate the possibility of a first-order phase transi¬ 
tion by solving Eq. (10) numerically. In Fig. (3), we show 
the non-equilibrium phase diagram in the T-oj plane for 
a fixed h/kjj. Starting from small drive frequencies, the 
line of second-order phase transitions reaches a tricritical 
point located at (w*(/i/w), T* = T® q ) and turns into a 
line of first-order transitions for larger w. 


Discussion. Besides the demonstration that driven- 
dissipative conditions can strongly reshape the phase di¬ 
agram of the Ising model, this study allows us to shine 
a new light on the fate of the universal properties of 
this model and, by extension, other similar models when 
driven to non-equilibrium steady states. When the drive 
is finite, we have found that the critical exponents (and 
the critical temperature) are strongly dependent on the 
details of the bath, thus loosing much of their universal¬ 
ity. We hope to report soon on the influence of dimen¬ 
sionality (i.e. away from the mean-field approach) on 
these results by studying the one-dimensional case via 
exact methods j351 . 
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Supplementary Material 
Stability of the mean-field solutions 


Here we check whether the non-zero mean-field solu¬ 


tions in Eq. (14) are stable. We start with a Master 
Equation for the probabilities P| and in terms of the 
rates P|f and P-|q: 

Pf Pf 4 Pf + P44- Pf 

Pf = + Pf, Pf — Pit Pf ■ 

Using Pj. = (l±y>)/2 and P^ = {\zpp)/2 for the ferromag¬ 
netic and anti-ferromagnetic cases, respectively, yields 

±ip = [Pf t M - PtiM] - [PitM + PfiM] (=M 

or, equivalently, 


t P = ~ [Pfl-M + PfiM] 1 <P -F 


P+tfc) - Pti(^) 

PitM + PtiM 

The quantity in curly brackets vanishes at the station¬ 


ary point, and gives precisely the condition in Eq. (10). 


Let p be this stationary point solution. To consider the 
stability of fluctuations, we expand p = Cp + Sip. The ex¬ 
pansion of the terms in curly brackets start at order Sip 
(because p is where it vanishes); so to lowest order, the 
term in square brackets does not need to be expanded. 
The linearized stability equation becomes 


where 


Sp = - [Ppt-(y>) + PnM] [1 T C{p)\ Sp , 


c(<3) = — ( Pit (<p)- R n( < p) \ 
dp \R lt (p) + R n (p) J 


Notice that R^(p) + R^i(p) > 0, so the stability of the 
solution rests upon whether [1 =F C(p)\ > 0. 

Using Eq. © , we find 


1 =F C(p) = =F /3zJ 


d K\2zJp\ s + A 


dp K\2zJp\ s 


B 


= =p fdzJ (A - B) 


1 


dp K\2zJp\ s 


B 
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The quantity in the square bracket above is al¬ 
ways negative. Therefore, the sign of 1 =F C((p) is 
that of ±sgn(J) (A — B). Now recall that T c = 
±sgn(J) [A/B(T C )\ T® q is larger than T® q for the 


non-trivial magnetization to be well defined; therefore 
±sgn( J) A > B. Thus, ±sgn( J) (A - B) > B \ :: 
sgn( J)] > 0. Hence, we conclude that the sign of 1=F C((p) 
is positive and the solutions we found are stable. 



